Angular-dependent upper critical field of overdoped Ba(Fei_ x Ni x ) 2 As 2 
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In-plane resistivity measurements as a function of temperature, magnetic field and its orientation 
with respect to the crystallographic ab— plane were used to study the upper critical field, if C 2, of 
two overdoped compositions of the iron-based superconductor Ba(Fei_ cc Ni cc )2As2, x =0.054 and 
x=0.072. Measurements were performed using precise alignment (with accuracy less than 0.1°) of 
magnetic field with respect to the Fe- As-plane. The dependence of the H C 2 on angle between the 
field and the afr-plane was measured in isothermal conditions in a broad temperature range. We 
found that the shape of H C 2(9), clearly deviates from Ginzburg-Landau functional form. 

PACS numbers: 74.70.Dd,72.15.-v,68.37.-d,61.05.cp 



I. INTRODUCTION 

The theory of the upper critical field in superconduc- 
tors was essentially developed by the mid-GOs, 1 3 with 
the linear H C 2(T) behavior close to T c described by 
the anisotropic Ginzburg-Landau theory leading in 
tetragonal crystals to dependence of H C 2 on angle with 
respect to the afr-plane: 
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Renewed interest to studies of the upper critical field 
was brought by the discovery of materials with poten- 
tially unconventional pairing mechanisms, heavy fermion, 
organic and especially cuprate superconductors, with 
the well documented d-symmetry of the superconduct- 
ing order parameter in the latter case. 6 It was quickly 
recognized that the angular dependence of the upper 
critical field, in particular on the angle <p in the basal 
plane of the tetragonal crystals, can be used to probe 
the anisotropy of the superconducting order parameter P 
These ideas were further developed for unconventional 
superconductors with various exotic order parameters 
in a series of papers by K. Maki and co-workers.™^ 
They stimulated a series of experimental studies, in par- 
ticular in S^RuO^j^HUl heavy fermiorP^ and organic 
super conduct or s.HZHSl 

Although these developments were correctly catching 
the importance of the anisotropy of the superconducting 
gap structure for the angular-dependent H C 2, the the- 
oretical data analysis was oversimplified by assumption 
of simple cylindrical or spherical Fermi surface shapes. 
The importance of the Fermi surface topology for the 
anisotropy of the H C 2 was brought to focus by the discov- 
ery of pronounced multi-band effects in superconductiv- 
ity of MgE>2. 22 24 The upper critical field of this layered 
material is anisotropic and can be varied significantly by 
carbon and aluminum dopingP^HUl controlling the mean 
free path of the carriers and changing the inter-band cou- 



pling. Theoretical modeling explicitly included these ef- 
fects into consideration and produced the angular depen- 
dence of the H C 2(0) which was virtually identical to the 
dependence of Eq. [IJ 28 

Recent developments in understanding of the upper 
critical field were greatly stimulated by the discovery 
in 2008 of the layered FeAs superconductors, 29 which 
opened up a new avenue in high T c research. The up- 
per critical fields of iron pnictides are very highp^ and 
besides the potential for high field applications this 
brings the possibility of the paramagnetic effects at low 
temperatures J 3 - 2 ^ 

Since Fe-As layers form a main building block of all 
iron-based superconductors pUEH these compounds show 
anisotropy of the electronic structure, reflected in the 
anisotropy of the upper critical field. Unlike sig- 
nature layered materials, organic super conduct or s^H and 
cupratesp^ the values of in iron arsedines are small, 
for most compounds in the range 2 to 4ESJHEI w itfi few 
except ions^MH we re values as high as 7 to 10 are found, 
see Ref. 52 for review. Furthermore 7^ in iron ar- 
senides shows notable temperature dependence, presum- 
ably reflecting multi-band effects! 4 ^ Low anisotropy val- 
ues suggest that at least close to the critical in zero-field, 
T C (H = 0), we are dealing with orbital H C 2. 

The detailed study of the anisotropy of the upper crit- 
ical field along principal a and c directions as a func- 
tion of doping was undertaken in the series of electron- 
doped compounds Ba(Fei_ a; Co :c )2As2 (BaCol22 in the 
following). 42 It was found that 7^ changes significantly 
between under-doped x < 0.08 and overdoped regions 
of the phase diagram, presumably reflecting change of 
the electronic structure. Comparison with optimally 
doped Ba(Fe i_ a; Ni a ;)2As2 (BaNil22 in the following), 
x=0.04d^ 4 ^ suggested that the anisotropy is not par- 
ticularly sensitive to the amount of disorder introduced 
by dopant x. 

On the other hand, several studies pointed out complex 
superconducting gap structure of iron arsenides, both in 
BaCol2#^IU an d analogous BaNil22P For both these 
materials, it was suggested that notable superconducting 
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FIG. 1. (Color Online) Temperature-dependent resistivity 
of two samples of BaFei_ x Ni x As2 used in this study, with 
x=0.054 (slightly overdoped) and x=0.072 (strongly over- 
doped), with doping level indicated with arrows with re- 
spect to temperature-doping phase diagram of BaNil22 after 
Ref.PTHand 48 shown in the inset. Note pronounced curvature 
of the p(T) for T > T c , typical of overdoped compositions.^ 
Sample resistivity value is defined with accuracy of about 20% 
due to uncertainty of geometric factors, see Refs. [46] and [62] 
for details. 



gap modulations along c-axis should be characteristic of 
the overdoped regime. Since this gap variation might be 
reflected in the angular dependence of the H C 2(9), we de- 
cided to perform systematic study of these compositions. 
To our knowledge, there was only one study addressing 
angular variation of H c2 in iron-based superconductors p3 
providing no high angular resolution data. 

In this paper we perform a detailed characterization 
of the upper critical field in overdoped BaNil22 as a 
function of temperature and field direction. Our main 
experimental finding is a clear deviation of the angular 
dependence from Eq. [I] particularly strong for the di- 
rections of magnetic field close to H || c, where orbital 
effects should be the strongest. 



II. EXPERIMENTAL 

Single crystals of Ba(Fei_ ;E Ni :c )2As2 were grown from 
FeAs flux using a starting load of metallic Ba, FeAs 
and NiAs, as described in detail elsewhere.ESl Crys- 
tals were thick platelets with large faces correspond- 
ing to the tetragonal (001) plane. Actual composition 
of the crystals was determined using WDS x-ray elec- 
tron probe microanalysis PHI The two compositions stud- 
ied were on the overdoped side of the phase diagram, 
slightly overdoped x=0.054 (T c =16 K) and strongly over- 
doped £=0.072 (T c =7.5 K), whereas maximum T c =19 K 
is achieved at optimal doping, x opt =0.046j 47 1 48 1 see dop- 
ing phase diagram in inset in Fig. [I] 
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FIG. 2. (Color Online) During experiments in single axis ro- 
tation system of 35 T magnet, the direction of magnetic field 
was aligned parallel to the conducting plane by resistivity 
measurements in field H slightly lower than H c2 \\, in which 
sample resistance shows strong angular dependence, black line 
in the top panel. The curve was measured in one-sided mo- 
tion of the rotator to avoid backlash, with deep minimum 
corresponding to if || ab or 0=0 condition. The red open sym- 
bols show alignment measurements, taken in a second angular 
sweep of the same rotation direction, and stopped at 0=0. H 
and T sweeps were used to determine the phase diagrams in 
H || ab condition, and then magnetic field angle with respect 
to the plane was changed by continuing rotation of the sam- 
ple in the same direction as during alignment. Because the 
orientation of the sample in the third direction, perpendicular 
to the rotation plane, was set by eye there may exist non-zero 
angle ip between the field-rotation plane and the plane of the 
normal to the sample. In most cases this angle should be less 
than 5°. 



Samples for in-plane resistivity, p, measurements were 
cleaved with a razor blade into rectangular strips with 
typical dimensions, 2 x (0.1 - 0.3) x (0.03 - 0.1) mm 3 
and the long side corresponding to tetragonal a-axis. All 
sample dimensions were measured with an accuracy of 
about 10%. Contacts to the samples were made by at- 
taching silver wires using ultrapure tin, resulting in an 
ultra low contact resistance (less than 10 /i^). 63 Resis- 
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tivity measurements were made using a standard four- 
probe technique, producing the p(T) curves as shown in 
Fig. [I] After initial preparation, samples were charac- 
terized in PPMS system, and then glued by GE-varnish 
to a plastic platform, fitting single axis rotator of the 
35 T DC magnet in National High Magnetic Field Lab- 
oratory in Tallahassee, Florida. Sample resistance was 
checked after mounting and found to be identical to the 
initial value. High-field measurements were made in He- 
cryostat with variable temperature control inset (VTI) 
allowing for temperatures down to 1.5 K. 

The stepping motor driven rotator enabled in situ ro- 
tation with 0.05° resolution around a horizontal axis in 
single axis rotation system of vertical 35 T magnetic field. 
During this rotation the direction of magnetic field with 
respect to the crystal stays in a plane of rotation, see 
Fig. |2j We can precisely align the direction of the mag- 
netic field parallel to the sample plane within the rotation 
plane, defined as 6 =0, using angular dependence of re- 
sistivity, measured in magnetic field slightly below H c2 \\. 
This alignment is illustrated in Fig. [2] In an ideal case of 
the second sample axis coinciding with the rotation axis, 
field-rotation plane should contain c-axis of the sample. 
There may have been a non-zero uncontrolled angle ip be- 
tween the field-rotation plane and the plane of the normal 
to the sample, see Fig. [2] We estimate that (p < 5°. 



III. RESULTS 

In Figs. [3] and [4] we show raw p(T) data for a set 
of magnetic fields aligned approximately along c-axis 
(0 =90°, top panels) and precisely along the conduct- 
ing plane (9 =0°, bottom panels), for BaNil22 samples 
with x=0.054 and x=0.072, respectively. We show also 
the lines corresponding to 20, 50 and 80% of resistiv- 
ity value immediately above the transition, p(T c ), used 
as criteria to determine the transition temperature as a 
function of magnetic field and construct the phase dia- 
grams, bottom panels (c) of Figs. [3] and [4] The use of 
these criteria is justified by small variation of the resistive 
transition width on application of magnetic field, and its 
independence on the extrapolation, typical problem for 
onset and offset criteria. 

As can be most clearly seen from the bottom panel of 
Fig. [3j the shapes of the H C 2 (T) phase diagrams in par- 
allel and perpendicular field orientations share the same 
features as found in previous studies of other Fe based 
systems. The H C 2, a b(T) flattens at low temperatures, 
while H C 2, C (T) maintains positive curvature down to the 
lowest temperatures of our experiment. Both these fea- 
tures are typical for layered materials, see for example 
Refs. 37 and M 

In Figs. [5] and [6] we show field dependences of in-plane 
resistivity taken at fixed temperatures with inclination 
angle as a parameter for slightly overdoped sample with 
£=0.054 and strongly overdoped sample £=0.072, respec- 
tively. The data analysis will be presented in the next 
section. 
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FIG. 3. (Color Online) In-plane resistivity p a vs. temper- 
ature for slightly overdoped Ba(Fei_ cc NL c )2As2, x=0.054 in 
magnetic fields (a) parallel to the conducting ab plane (right 
to left, T, 1 T, 2.5 T, 5 T, 7.5 T, 10 T, 15 T, 17.5 T, 
20 T, 25 T, 30 T, 35 T); (b) parallel to the c-axis (right to 
left, T, 1 T, 2 T, 4 T, 8 T, 10 T, 12 T, 15 T, 20 T, 30 T, 
35 T). Lines indicate 20, 50, and 80 % of the resistivity value 
immediately above the superconducting transition. Bottom 
panel (c) shows H c2 (T) phase diagrams for both directions of 
magnetic field. 



IV. DISCUSSION 
A. Angular dependence of H C 2 

1. Linearization of the H C 2{0) dependence 

To check if Eq. [I] describes our data, instead of com- 
monly used data fitting, as shown in the bottom panels 
of Figs. [5] and [6j we used an approach based on data 
transformation so as to make possible deviations clearly 
visible. According to Eq.[l] the H~ 2 2 vs (sin 2 0) should be 
a straight line, and in Figs. [7] and [8] we plot this way the 
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FIG. 4. (Color Online) In-plane resistivity p a vs. temper- 
ature for heavily overdoped Ba(Fei_ cc Ni cc ) 2 As2, x=0.072 in 
magnetic fields (a) parallel to the conducting ab plane (right 
to left, T, 1 T, 2.5 T, 5 T, 7.5 T, 10 T, 15 T, 17.5 T, 20 T, 
22.5 T); (b) parallel to the c-axis (right to left, T, 1 T, 2 T, 
4 T, 6 T). Lines indicate 20, 50, and 80 % of the resistiv- 
ity value immediately above the superconducting transition. 
Bottom panel (c) shows H C 2(T) phase diagrams for both di- 
rections of magnetic field. 



data for samples with x =0.054 and x =0.72 respectively. 
The data show clear deviation from linear trend, irre- 
spectively of the criterion of H C 2 determination from the 
resistivity data, with the deviation being the strongest 
close to H || c— axis or sin 2 0=1. To check if the devia- 
tion from Eq. [I] in Figs. [7] and [8] can be caused by finite 
inclination angle cp (see Fig. [2] for the definition), here we 
provide the angular dependence of H C 2 for arbitrary ip 
Choosing the cross-section of the plane, in which H is ro- 
tated, with the ab crystal plane (see Fig.|2| as the x axis, 
we obtain in the crystal frame c = (0, 0,1) and the unit 
vector along the field h = (cos#, sin#sin(/?, sin#cos(/?). 
This gives for the angle C between the field and c axis: 
cos# c = c ■ h = sin cos cp. We then obtain for geometry 



FIG. 5. (Color Online) Field dependence of in-plane resistiv- 
ity p(H) of slightly overdoped Ba(Fei_ a; Ni ;r )2As2, x=0.054 
sample at T=13 K (panel (a)) and T=9 K (panel (b)) 
with magnetic field inclination angle as a parameter, (c) 
Isotherms H c2 (0), obtained at 9 K and 13 K, using 80%, 50% 
and 20% criteria. Solid line shows fit to Eq. [I] 



of our experiment: 
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It is seen that constant (/?, as determined by our ex- 
perimental geometry, does not change the linear relation 
of H~2 vs sin 2 despite changing the magnitude of the 
variation, vanishing for cp = 90°, corresponding to field 
rotation parallel to the conducting plane. Therefore the 
linear dependence of on sin 2 is not affected by a 
misalignment (p. 

The H c2 (0) described by the Eq. [TJ is a direct con- 
sequence of the linearized GL equation for anisotropic 
materials at H C 2'. 



(e 2 )»fcii < n fc * = *, 



(3) 



5 



x=0.072 




0.00 0.25 0.50 0.75 1.00 



0.008 



H\\ab 




30 60 
(degrees) 



FIG. 6. (Color Online) Field dependence of in-plane resistiv- 
ity p(H) of strongly overdoped Ba(Fei_ :c NL c )2As2, x=0.072 
sample at T =1.5 K (panel (a)) and T=4 K (panel (b)) 
with magnetic field inclination angle 6 as a parameter, (c) 
Isotherms H c2 (0), obtained 1.5 K and 4 K, using 80%, 50% 
and 20% criteria. Solid line shows fit to Eq. [I] 



where II = V+27riA/(/>o, A is the vector potential and 0o 
is the flux quantum; summation is implied over repeating 
indices. Both sides of this equation are scalars, so that 
(£, 2 )ik is a second rank tensor with the standard angular 
dependence which is reflected in Eq. ([I]). 

We note that in the original papers, the angular de- 
pendence, Eq.Q, has been derived for single band s- 
wave superconductors. It has also been recently shown 
that this behavior is expected for arbitrary Fermi sur- 
face, the superconducting gap modulation and for multi 
- band materials.^ However, this conclusion is achieved 
assuming the explicit factorization of the pairing poten- 
tial and order parameter, V(fe, fe') = VoQ(k)Q(k') and 
A(T, k) — \£(r, T) f2(kp), see detailed discussion in the 
theory section below. There is no microscopic justifi- 
cation for such factorization in complex superconductors 
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FIG. 7. (Color Online) Angular dependence H C 2(9), deter- 
mined from fixed temperature p(H) of Fig. [5] using 20%, 
50% and 80% criteria (top to bottom), for slightly overdoped 
Ba(Fei_ cc NL c ) 2 As2, x=0.054 at 9 K (top panel) and 13 K 
(bottom panel). The lines are guides to the eyes. The data 
are plotted as if ^ 2 (sin 2 0), which according to Eq. [I] should 
be a straight line. The lines are guides to the eyes. 



and deviations from Eq.Q can be naturally explained by 
violation of this procedure. In addition, for iron-pnictides 
the importance of the paramagnetic effects for magnetic 
fields parallel to the Fe-As pl ane was suggested to explain 
unusual shape of the H c2 (TW^^). This may also lead 
to the deviation from Eq.dTl) with the maximum effect 
expected at low temperatures and for orientations close 
to H || ab— planes. 

Motivated by these considerations, we compile in Fig. [9] 
the published data for various layered materials, ana- 
lyzed by plotting H~ 2 2 vs sin 2 0. The data are arranged 
with decreasing anisotropy from top to bottom. The 
most anisotropic materials, staged graphite intercalation 
compounds (top panel, data from Ref. [66j) and layered 
Sr2Ru04 (data from Ref. 15) closely follow Eq. [I] In- 
terestingly, that clear deviations from this behavior in 
Sr2Ru04, arising due to unusual limiting mechanism 
in magnetic fields close to H \\ ab^ is very difficult 
to recognize in a limited angular range near #=0, as 
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FIG. 8. (Color Online) Angular dependence H C 2(9), deter- 
mined from fixed temperature p(H) of Fig. fusing 20%, 50% 
and 80% criteria (top to bottom), for strongly overdoped 
Ba(Fei_ cc NL c ) 2 As2, x =0.072 at 1.5 K (top panel) and 4 K 
(bottom panel). The lines are guides to the eyes. The data 
are plotted as if ^ 2 (sin 2 0), which according to Eq. [I] should 
be a straight line. 



the dependence in the whole range is dominated by the 
anisotropy of the Fermi surface. On the other hand, 
two materials in which superconductivity shows strong 
multiband features, MgBj^and NbSe2 67 show distinctly 
different angular dependences. The H~^(sin0) in pure 
MgB 2 68 shows downward bent as field approaches c— 
axis, #=90°, similar but much less pronounced than in 
our observations in BaNil22. On the other hand, doped 
Mg(Bi_ x Al a: )2 closely follows the linear i^ 2 (sin#) de- 
pendence, Eq. [T] which may be suggestive that doping 
diminishes multi-band effects due to enhanced inter-band 
scattering. For pure NbSe2 the i^ 2 (sin#) plot shows 
most clear deviations from linearity among all materials, 
with an upward curvature towards 6 = 90°, of oppo- 
site trend to pure MgB2 and BaNil22. The two angular 
data sets for profoundly multi-band iron pnictide super- 
conductors, slightly underdoped, BaK122, 60 and heav- 
ily overdoped K122, 65 generally follow linear dependence 
despite profound difference in the superconducting gap 



structure, nodeless in the former case^and with vertical 
line nodes in the latter .^Considering that among all the 
materials for which we were able to find published H C 2 (0) , 
only pure multi-band MgB2 and NbSe2 reveal clear devi- 
ations from Eq. 0, it is tempting to relate the observed 
deviations with the multi-band superconductivity in the 
clean limit. This might be quite natural that in these 
systems the factorization of the pairing potential and of 
the order parameter does not hold given the complexity 
of the in- and inter - band interactions. This explana- 
tion, however, is not universal, since multi-band effects 
are very pronounced in high purity crystals of KFe2As2, 
but no clear deviations from Eq.Q are found there. On 
the other hand, it is hard to consider overdoped BaNil22 
as a clean system, since scattering due to substitutional 
disorder, especially on the Fe site, is significant in these 
compositions. The observation that the deviations from 
the linear plot in MgB2 diminish with disorder, suggest 
that it is /c-dependence of the gap magnitude, rather than 
multi-band nature of the Fermi surface itself is important 
for the unusual angular dependence. This conclusion is in 
line with the recent extension of the HW theory for multi- 
band superconductors with arbitrary Fermi surfaces.^ 

In discussing these results we should keep in mind, that 
in all cases, except for Sr2Ru04, the H C 2 was measured 
resistively, so that inevitably its determination is approx- 
imate since the resistive transition as a rule has finite 
width and hence, the H c2 values depend on a criterion 
chosen. Finite resistivity in the flux - flow regime (most 
pronounced in the clean systems) broadens the transition 
making resistive determination difficult. From this point 
of view, assertions of Kim et alP^ that their data allow 
to distinguish between two models, GL and two - band 
Usadel approach by Gurevichf^ are hard to accept. 

In compounds with relatively high T c , the determina- 
tion of H C 2 from resistive measurements is also compli- 
cated by the phenomenon of vortex lattice melting: above 
the melting point, the resistivity is close to that of the 
normal phase and H C 2 per se becomes invisible in resistiv- 
ity measurements. This complication in a given material 
might affect the measurements stronger near T c than at 
low temperatures. 



2. Theoretical considerations 

Recently, the theoretical justification for the near- 
universality of the angular dependence (Eq. [I]) has been 
offered in Ref. 71 for orbital H C 2{9) albeit under certain 
restricting assumptions. We reproduce briefly the argu- 
ment of this work. Some deviations from the effective 
mass model have been predicted theoretically for quasi- 
one-dimensional conductors , see Ref. [72] We also note 
that the theory of H C 2{6) at the arbitrary angle with re- 
spect to the crystal axes in multi-band superconductors 
has additional complications. In particular, if H is not 
parallel to the c— axis and the mass anisotropy ratios in 
two bands are different, and the first Landau level solu- 
tions which are then substituted into the HW integrals 
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cannot simultaneously satisfy the two - band Eilenberger 
equations and the gap eq uatio ns, and the higher Landau 
levels become important ! 52 * 76 * Perhaps this provides an 
alternative explanation of why we observe the deviations 
from the conventional scaling only in clean two - band 
systems.^ Keeping this in mind, we proceed with the 
generalized description! 7 -^ 

We recall that, as HW^ had shown within the weak- 
coupling BCS-Gor'kov theory, the linearized GL equation 
— £ 2 II 2 A = A for the isotropic case holds at any T at 
H c2 (T). However f (T) does not diverge for T T C (H) < 
T c (0), -instead it is finite and should be evaluated to 
satisfy the basic self-consistency (or gap) equation for the 
order parameter A. It is easily shown then that H C 2(T) = 

For materials with anisotropic Fermi surfaces and or- 
der parameters one can consider Eq. ([3| as an ansatz for 
satisfying the self-consistency equation of the theory. To 
facilitate this derivation one assumes the coupling poten- 
tial of the form V(fe, k') = VoQ(k)Q(k') and the order pa- 
rameter A = Sfr(r,T) Q(k F ). The function Q(k F ) which 
determines the kp dependence of A is normalized so that 
(f} 2 ) = 1. This convenient approximation has been first 
designed by Markovic and Kadanov^ and it works well 
for one band materials with anisotropic coupling. 

For two-band materials one can use a generalization 
of this procedure. Other than s-wave order parameters 
emerge if the coupling V(fe, k') responsible for supercon- 
ductivity is not a constant (or a 2 x 2 matrix of k indepen- 
dent constants for two bands). The formally simplest way 
to consider different from s-wave order parameters with- 
out going to details of microscopic interactions is to use 
a "separable" potential, V a p(k, k') = V^ft a (k)ft^(k f ) 

where is a k independent matrix, and look for the or- 
der parameters in the form A a (r , T, fe) = (r , T)Vt a (k) 
with the normalization (Q 2 ) a = 1 for each band (a, /3 = 
1, 2 are band indices). 

It has been shown in Ref. [7T]that in this approximation 
the ansatz 

- (fWiKkka = A a , a = 1,2. (4) 

leads to solutions satisfying the self-consistency equa- 
tions, which provide all components of the tensor 
(£ > 2 )ik(T). In this way the ansatz Q is proven cor- 
rect. Along with (£ 2 )^(T), the upper critical field and 
its anisotropy can be evaluated for a given Fermi surface 
and the order parameter symmetries Q a . 

Since (£, 2 )ik is a second rank tensor, we conclude, ex- 
actly as in the case of anisotropic GL Eq. (J3|, that the an- 
gular dependence of Eq. should hold at any tempera- 
ture and for any anisotropic Fermi surface, any symmetry 
of the bands order parameters, and at any temperature. 
In other words, within the approximation of "separable" 
coupling potential, the dependence ([!]) is "universal" . 

Clearly, "separable" potentials do not exhaust all pos- 
sible interactions and, therefore, other forms of the angu- 
lar dependence H C 2{6) can exist. An example of such a 
potential has been studied in Ref. [9] Such potentials may 



lead to gradient terms in GL equations different from the 
standard forms (|3| or Q and, therefore, different from 
Eq. ([I]) angular dependencies, see e.g. Ref. [751 

We therefore may conclude that deviations of the ob- 
served angular dependence of H c2 from the form ([!]) (or 
deviations of H~ 2 2 plotted vs. sin 2 from the straight 
line) signal that the coupling potential cannot be written 
in the separable form. On the other hand, the exam- 
ple of separable potentials (for any Fermi surface and 
any order parameter symmetry) shows that there is no 
direct relation between the angular dependence of H c2 , 
Fermi surfaces, and order parameter symmetries. How- 
ever, deviations of H C 2{6) from the form may carry 
such an information. To investigate this question fur- 
ther one would need a better data on these deviations, in 
particular, criterion-independent determination of H c2 , 
which is hard to achieve in resistive measurements. On 
the theoretical side, of course, one should go beyond the 
weak coupling and separable coupling potentials. 



B. Temperature dependence of H C 2 

There are two mechanisms that determine the up- 
per critical field of superconductors. The first one, de- 
termined by the supercurrent flow to screen the mag- 
netic field, is referred to as orbital limiting and de- 
scribed by HW theory The upper critical field at 
T — >• limit, if C 2(0), in HW theory is determined by 
the slope of the H c2 (T) curve close to T c , and as T goes 
to zero the curve shows downward deviation from linear 
dependence and eventual saturation towards the value 
H c2 (0) « 0.7T c ^f in isotropic case. 

Rather rare exceptions, when the upper critical field is 
not determined by the orbital motion, are found in the 
materials in which orbital motion of electrons is ham- 
pered by either short mean free path, heavy mass of 
conduction electrons in heavy fermion materials or weak 
links between the conducting layers in Josephson struc- 
tures or in naturally layered materials, provided that the 
magnetic field is aligned precisely parallel to the con- 
ducting layer. In this situation the upper critical field 
H c2 is determined by Zeeman splitting of electron levels, 
known as Clogston-Chandrasekhar 3 paramagnetic limit. 
This field is determined by a decrease of paramagnetic 
energy becoming equal to condensation energy of super- 
conductor. In weak coupling BCS superconductors the 
paramagnetic limiting field is determined in T — »• limit 
as H P =1.S4T C , where H p is field in Tesla and T c is in K. 

As can be seen from Figs. [3] and [4] the upper criti- 
cal fields in H || ab configuration are higher than the 
weak-limit paramagnetic limiting H p , equal to 32.2 T 
(x =0.054) and 13.8 T (x =0.072). These high values 
may come from strong coupling nature of superconduc- 
tivity in iron pnictides, or indeed reflect paramagnetic 

limitin g at low temperatures as was suggested in several 
studiesPMEU 
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CONCLUSIONS 



anisotropic pe netratio n depth and thermal conductivity 
measurements pfl 54 l 59 l 



By performing high angular resolution study of the 
upper critical field in two overdoped compositions of 
iron pnictide superconductor Ba(Fei_ iC NI E )2As2, we find 
clear deviations from the anisotropic Ginzburg-Landau 
form. Implementing linearization plot analysis of our 
and previously published data, we find clear deviations 
from the form only in the case of clean fully - gapped 
multi-band superconductors, such as NbSe2 ang MgB2, 
but not in dirty MgB2 and not clean KFe2As2- We 
speculate, that the dependence may reflect c-axis mod- 
ulation of the superconducting gap, as suggested by 
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FIG. 9. (Color Online) Analysis of the isothermal angular dependence of H C 2 on inclination angle to the highly conducting 
plane using linearization plot if ^ (sin 2 9). Left panels show digitized H C 2(0), right panels plot the same data as if J^ 2 (sin 2 0): 
(a) Graphite intercalation compounds 66 C4RbHg (T c =0.99 K, measurements taken at Xh=0.44 K, open circles) and C4KHg 
(T c =0.73 K, T h =0.40 K, solid squares); (b) Sr 2 Ru0 4 (T c =1.43 K, T^=0.10 K, Ref. 15); (c) Mg(Bi_,AU)2, Ref . EH x=0.12 
(T c =30.8 K, T h =U K, black solid squares, and T h =23 K, red solid circles) and x=0.21 (T c =25.5 K, TWO K, blue open 
circles); (d) NbSe 2 , Ref. EH (T c = 7.2 K, T h =4.2 K); (e) (Bai-^K^FesAss, Ref. 60, (T c =28 K, T h =20 K, using different criteria 
for resistive transition, zero resistance- black triangles, midpoint- red circles, onset - green squares); (f) KFe2As2, Ref. [65] 
(T c =3.8 K, 7V=0.5 K). 



